Can LSND be included in a 3-Neutrino framework? 
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We study the special features emerging from a three lepton generation 

analysis of the available neutrino oscillation data (solar, atmospheric and 

LSND). First we find that it is possible to explain all three sets of data in terms 

of the standard left handed neutrinos without the need of sterile neutrinos. 

Second we find a significant difference in the mass matrix extracted from the 

data, depending on the analysis (without or with LSND), if the mass of the 

lightest neutrino, which cannot be determined from the neutrino oscillation 

data alone, is relatively small, i.e <0.1 eV. To compare with other processes we 

used the R-parity violating Minimal Supersymmetric Standard Model (^? p - 

MSSM) for the theoretical description of the neutrino masses. Using the 

oscillation data we were able to constrain the parameters of the model. In 

particular we were able to obtain values for the coupling constants of the 

$ P -MSSM. 

PACS number (s):{} 

Searches for neutrino masses have intensively been performed during the last decades. 
But till now no evidence for massive neutrinos was found in experiments aiming to measure 
the mass directly, but only upper limits were set as follows 

m Ve < 3.5eV,m^ < 160KeV, m„ T < 23MeV. (1) 

On the other hand some hints for massive neutrinos have been seen more than 30 years 
ago in experiments measuring the electron neutrino flux coming from the sun ||. In these 
experiments the flux of u e coming from the sun to the earth was found to be much smaller 
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than the expectations. This observed difference of the expected vs the measured flux could 
be explained in terms of neutrino oscillations. The v e mixing, due to oscillations to other 
neutrino flavors, leads to a smaller number of u e 's at the detector. These oscillations between 
different states can only occur if the weak eigenstates are different from the mass eigenstates 
and if the mass eigenstates are not all degenerate. This means that some neutrinos must be 
massive. The probability for a v, which was produced in the flavor state \v a ) with energy 
E, to be detected in the flavor state \vp) is given in a three family mixing scheme by 



P(a ->/?) = 8 a ,p - 4 u *iUpiU aj Up.j sin 2 
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where L is the source detector distance. Here Am?- = \mf — m 2 | is the difference of the 
squared masses of the mass eigenstates i and j, U a i are the elements ai of the mixing matrix 
U, assuming it is real (negligible CP violation), which connects the flavor eigenstates \v a ) 
and the mass eigenstates by |z/ Q ) = X^KnK)- If CP- violation is negligible the mixing 
matrix U can be parameterized in analogy to the CKM matrix by three angles, 

^ C12C13 S 12 Ci3 s 13 ^ 
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-S12C23 - C12S23S13 C12C23 - S12S23S13 S23C13 
S12S23 _ C12C23S13 -C12S23 - S12C23S13 C23C13 

where and Cjj stand for sin(%) and cos(%), respectively. 

The search for neutrino oscillations has been the subject of many recent experiments 
PHlOU- The experiments which found evidence for the existence of v oscillations can be 



characterized by the sources of neutrinos which they use as solar atmospheric |j-|8| and 
accelerator neutrino experiments |9|,|lIJ] . For the atmospheric and solar neutrino experiments 
several groups find evidence for neutrino oscillations. In the field of accelerator neutrino 
experiments only one group, the LSND collaboration, has claimed evidence for neutrino 
oscillations, but this result has not yet been confirmed by any other experiment. Furthermore 
it is commonly believed that all the three experiments (solar, atmospheric and LSND) cannot 
be accommodated with just two independent Am 2 . Thus the LSND results were put in 



doubt. We will see that this outcome is due to the essentially two flavor scenario, which has 
hitherto been employed in most analyses. We will see that in a rigorous three generation 
analysis all the available data can be accommodated and can be used to constrain the 
neutrino mass matrix. Since the emerging solution is, unfortunately, not unique, we propose 
a way to test in the future this solution and thus, indirectly, to test the LSND results. 

As we have already mentioned, the analysis of neutrino oscillation experiments is of- 
ten done independently for each experiment within a two family scenario. The oscillation 
probability with the same notation as in eqn. @ in a two generation mixing scheme is 

P( a -+p) = sin 2 26 sin 2 ^1.27Am 2 ^ , (4) 

where Am 2 is in Units eV 2 , L in m and E in MeV. The result of analyses of the solar neutrino 
experiments for the mass squared difference Am 2 by including matter effects (Mikheyev- 
Smirnov-Wolfenstein [pT|,[12| ) is 

4 x 10- 6 eV 2 < Am 2 sun < 1.2 x 10~ 5 eV 2 . (5) 

For the atmospheric v experiments one finds disappearance oscillations of probably into 
v T and obtains for Am 2 

4 x 10" 4 e1/ 2 < Am 2 < 6 x 10~ 3 eV 2 . (6) 

Fitting now the results from the LSND collaboration for oscillations of — > 9 e coming form 
yU + decay at rest and — > v e which originate form the decay of 7r + also in a two family 
analysis yields a mass splitting of 

O.lel/ 2 < Am 2 < l.Oel/ 2 . (7) 

The three different splittings Am 2 of eqns. (|) to (^) seem to be not compatible with a 
mixing of only three v families and led to the conclusion that either the LSND result is 
wrong or that there should be a fourth sterile neutrino which mixes with the others. But 
before one makes such severe claims one should carefully scrutinize the methods which have 
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been used to obtain these results. By looking at eqns. (|2]) and (|J) one realizes that the 
squared mass differences obtained by fits using a two family scenario and by using a three 
family scenario need not to be the same. Indeed the squared mass difference extracted 
form eqn. (|j) are a convolution of the squared mass differences in the three family scheme. 
Therefore the need for a sterile neutrino does not necessarily follow from the above argument. 
Thus one should analyze the experimental data not in the oversimplified two family but in 
a three family mixing scheme. This has been done by several authors [p^-|l9fl. The results 
of these analyses are displayed in table | and |T|. In table | we show the results of such 
analyses which excluded the LSND results while in table || the results of analyses including 
the LSND results are given. At first we want to examine the influence of the inclusion of 
the LSND result on the oscillation parameters in the different analyses. By comparing the 
two tables we see that the influence of the inclusion of the LSND experiment on the mixing 
angles #13 and #23 is negligible and is still relatively small for 812- The mass splitting Amf 2 of 
the lower lying masses is found to be larger by including the LSND result in most analyses. 
Also for the mass splitting Am^ the inclusion of the LSND result leads to a larger splitting. 
In average, the inclusion of the LSND result in global fits to all data leads only to a larger 
mass splitting Amfj, which is due to the small oscillation length (L ~ 20m) in LSND. To see 
now the consequences of this difference we calculated the entries of the mass matrix in the 
weak basis for the different analyses. Oscillations fix the difference of the masses squared 
but not the absolute scale of the masses. For values of the smallest mass m\ larger than 1 
eV we find no significant difference between the two types of analyses. For three cases of 
the smallest mass, m x = 0.0, 0.01 and 0.1 eV, we show the average value for the neutrino 
mass matrix for the two types of analyses in table |TTT| . By comparing the results one finds 
that there exists a significant difference between the two types of solutions for these small 
masses. Because of the non observation of the neutrinoless double beta decay (0^/3/5) one 



expects such small masses for Majorana neutrinos pO| , pT| . As it is well known the mass 
extracted from the 0z/ double beta decay experiments is the absolute value of the quantity: 
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< m v > 



E^Afm, (8) 



where A^ p are the CP eigenvalues of the neutrino mass eigenstate Vj with mass rrij, which, 
of course, is positive. From the data of TABLE II we find: 

< m v >= 0.006,0.013,0.101, respectively (without LSND) (9) 

< m u >= 0.026,0.032,0.115, respectively (with LSND) (10) 

The small values can differentiate between the two types of analysis, but, unfortunately, 
they are far below the present experimental limits. They are, however, within the goals of 
the planned experiments p0|^l| . The experimentally most interesting values obtained for 
large m\ do not seem to depend on the type of the analysis (with or without LSND). 

To decide now which of the two kind of analyses is correct one has to check the consistency 
of the results with other physical observables. To compare with other processes we have to 
consider now a model which is able to accommodate massive neutrinos. This means that 
we have to use a model which is an extension of the Standard Model (SM). One natural 
way to extend the SM is supersymmetry (SUSY) ||22j| . A popular SUSY extension is the 
R-parity violating supersymmetric SM (]?L p -MSSM) . Within this model we have for each 
SM particle a supersymmetric partner, called sparticles. For details about the ft p -MSSM 
see e.g. |[23|| . Within this model all neutrinos acquire masses by mixing with the the SUSY 
partners of the photon, the Z° and the two neutral SUSY Higgs fields, see fig. [[], and by 
loop corrections, see fig. 0. For details how neutrinos acquire masses within this model see 



we 



26] . Using the values for the entries of the neutrino mass matrix given in table [L I 
extracted the values of the coupling constants for single term dominance. They are given in 
table [TV]. We see that for the couplings exists a significant difference only in the cases 
of mi = 0.00 and m\ = 0.01 eV. These small masses are found in a predictive model which 
use an underling symmetry to describe quark and lepton masses |24[ . 



To conclude, we discussed the analyses of neutrino oscillation data and showed that it is 
necessary to use a three family mixing scheme to do such an analysis. We presented results 



of different three family analyses including and excluding the LSND result and found that 
the inclusion mainly influences the mass splitting Am| 3 . We found that to distinguish the 
two different approaches and test by this the LSND result will only be possible if the smallest 
mass is smaller than 0.01 eV. 
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Am 2 2 [eV 2 } Am 2 3 [eV 2 } 9 12 [Deg.] B 23 [Deg.} O^Deg.] ref. 

3 x 10~ 6 - 7 x 10" 5 0.01 53-62 28-37 <13 

4 x 10~ 6 - 7 x 10~ 5 1.0 51-72 27-32 <4 || 



0i 2 [Deg.} 


e 23 [Deg.] 


0is[Deg.} 


53-62 


28-37 


<13 


51-72 


27-32 


<4 


51-72 


28-33 


<3 


39.23 


45 


26.6 



4 x 10~ 6 - 7 x 10~ 5 0.1 51-72 28-33 <3 H 



10~ 4 8 x 10" 4 39.23 45 26.6 @ 



TABLE I. Results of neutrino oscillation analyses in a three family mixing scheme excluding the 
LSND results. 
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Am 2 12 [eV 2 } 


AmUeV 2 } 


6 l2 [Deg] 


#23 [Deg.] 


0n[Deg] 


ref. 


10" 4 - 10" 3 


0.3 


35.5 


27.3 


13.1 


1 


1§ 


io~ 4 - i(r 3 


0.3 


54.5 


27.3 


13.1 




15| 


2.87 x 10~ 4 


1.11 


45 


28.9 


4.2 




17| 


10" 4 - 10" 3 


0.4 


37.6 


26.5 


10.3 




16| 


4 x 10" 6 - 7 x 10~ 5 


1 


51-72 


27-32 


3-4 




13| 



TABLE II. Results of neutrino oscillation analyses in three family mixing scheme including the 
LSND results. 



mi = 0.00 eV 



mi = 0.01 eV 



m 1 = 0.10 eV 



^ .101 .005 .005 ^ 



.005 .181 .076 
.005 .076 .186 



without LSND 



.006 .012 .007 
.012 .109 .099 
.007 .099 .111 



eV 



.013 .010 .008 
.010 .114 .095 
.008 .095 .117 



eV 



eV 



LSND included 



.026 .064 .072 
.064 .292 .339 
.072 .339 .457 



eV 



.032 .061 .073 
.061 .296 .336 
.073 .336 .459 



eV 



.115 .050 .067 
.050 .354 .301 
.067 .301 .496 



eV 



TABLE III. Results for the neutrino mass matrix in the weak basis. Given are the averaged 
values for matrix elements using the analyses discussed in tables Q and |Il| for Xj =1 which we 
found to give the largest differences. With LSND one obtains larger off-diagonal elements. 



S 









without LSND 






LSND included 




mi [eV] 




0.00 


0.01 


0.10 


0.00 


0.01 


0.10 


v 

133 


x 10 4 


0.4 


0.6 


1.6 


0.8 


0.9 


1.7 


^M SUSY /100GeV 


A233 


x 10 4 


1.6 


1.7 


2.1 


2.7 


2.7 


2.9 


^Msusr/WOGeV 


\i 

A 333 


x 10 4 


1.6 


1.7 


2.1 


3.3 


3.3 


3.5 


y/Msusr/WOGeV 


-^133 


x 10 4 


1.7 


2.6 


7.1 


3.6 


4.0 


7.6 


y/M SUSY /100GeV 


A233 


x 10 4 


7.4 


7.5 


9.5 


12.1 


12.2 


13.3 


y/Msusv/WOGeV 



TABLE IV. Values for trilinear couplings for single term dominance using neutrino oscillation 
analyses and different mass schemes for mi. For the calculation of these values of the trilinear 
coupling constants we assumed as usual that masses and soft parameters are of the same order and 
are parameterized by Msusy which is expected to lie in the region from 100 to 1000 GeV. 
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FIG. 1. Mixing of the neutrinos and neutralinos. 



FIG. 2. Radiative contributions to the neutrino mass in the ^ P -MSSM. 
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Fig. 2 
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